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General  Summary 

^  Geomagnetic  storms  arise  from  a  large-scale  electrical  discharge  process 
surrounding  the  earth  and  It  has  been  known  that  they  are  one  of  the  most 
serious  natural  hazards  to  the  Defense  Systems .  One  possible  solution  to  this 
hazard  is:  (1)  to  forecast  accurately  the  occurrence  and  intensity  of 
geomagnetic  storms,  (11)  to  alert  all  warning  systems  which  would  fall  during 
geomagnetic  storms,  and  (ill)  to  mobilize  emergency  systems  to  supplement 
them.  Thus,  the  accuracy  of  forecasting  the  occurrence  and  Intensity  of 
geomagnetic  storms  is  of  utmost  importance.  Reliable  forecasting  would  be 
able  to  reduce  drastically  expenses  Involved  in  mobilizing  these  emergency 
systems. 

We  have  recently  made  Important  progress  in  the  forecasting  by  the 
discovery  of  a  particular  solar  wind  parameter  which  controls  the  development 
of  geomagnetic  storms.  This  quantity  e  is  given  by 
e  -  VB2  sin4  (|)  *q2 

where 

V  -  the  solar  wind  speed 

B  »  the  solar  wind  magnetic  field  magnitude 
9  »  the  polar  angle  of  the  solar  wind  magnetic  field  vector 
■  constant  (*■  7  earth  radii) 

Having  this  quantity,  our  tentative  criterion  for  the  defense  systems  which 

are  vulnerable  to  geomagnetic  activity  may  be  given  as  follows: 
e  ^  ^q18  erg/sec  -  normal  operation 
1 8 

e  ~  5  x  10  erg/sec  -  require  alert 

19 

e  ~  10  erg/sec  -  prepare  for  system  failure 

20 

e  >  10  erg/sec  -  system  failure 


Major  Results 

I.  Establishment  of  the  outline  of  the  geomagnetic  storm  prediction  scheme 
for  the  defense  systems  which  are  vulnerable  to  geomagnetic  activity 


Summary 

Several  recent  major  advances  in  solar-terrestrial  physics  have  made  the 
prediction  of  geomagnetic  storms  a  definite  possibility.  The  operation  of 
defense  systems  which  are  vulnerable  to  geomagnetic  activity  can  be  greatly 
strengthened  by  a  computer  system  which  displays  the  predicted  problem  areas 
on  a  computer  screen  on  a  continuous  basis  with  a  6-12  hour  lead  time.  We 
have  most  of  the  major  computer  codes  for  our  purpose,  although  some  have  to 
be  Improved  to  meet  our  specifications.  Thus,  one  of  our  main  efforts  should 
be  to  assemble  the  codes  for  our  particular  purpose. 

1.  Introduction 


It  has  become  apparent  that  the  implementation  of  a  geomagnetic  storm 
prediction  scheme  will  significantly  improve  the  capability  of  defense  systems 
which  are  vulnerable  to  geomagnetic  activity.  Specifically,  it  is  Important 
to  predict  the  location  of  the  equatorward  boundary  of  the  auroral  oval  and 
characteristics  of  the  ionospheric  trough  (such  as  the  depth  and  width)  with  _a 
6-12  hour  lead  time  on  a  continuous  basis  and  display  this  information  using  a 
dynamic  computer  display  system.  Such  a  prediction  scheme  has  become  a 
definite  possibility  by  several  recent  Important  advances  in  solar-terrestrial 
physics.  Including  the  finding  of  the  solar  wind-magnetosphere  energy  coupling 
function  e,  the  success  of  the  first  three-dimensional  simulation  of  disturbed 
solar  wind  conditions  caused  by  solar  flares  and  the  development  of  various 
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computer  codes  in  modeling  the  Ionosphere.  In  Sections  IV  and  V,  some  of  the 
details  of  the  recent  progress  are  reported. 

Therefore,  it  is  important  to  set  up  a  practical  numerical  prediction 
scheme  specifically  for  the  defense  systems  which  are  vulnerable  to  geomag¬ 
netic  activity.  The  scheme  consists  of  a  number  of  powerful  computer  codes. 
Inputs  for  the  computer  codes  are  various  solar  observations  and  earthbound 
satellite  observations  (as  will  be  described  below  in  detail),  and  the  outputs 
are  the  equatorward  boundary  of  the  auroral  oval  and  the  trough  conditions 
which  are  displayed  on  a  dynamic  display  system.  What  is  Important  at  this 
point  of  the  development  of  the  storm  prediction  scheme  is:  (i)  to  identify 
all  necessary  computer  codes,  (11)  to  develop  missing  codes;  (ill)  to  assemble 
all  the  codes,  (lv)  to  test  and  improve  the  scheme. 

2.  General  Prediction  Scheme;  Figure  1 
(1)  Input 

The  input  consists  of:  (i)  coronal  hole  monitoring  and  (11)  solar  flare 
monitoring. 

(1)  Coronal  hole  monitoring 

The  coronal  hole  can  be  monitored  by 

(a)  The  white  light  coronal  observation  at  the  eastern  limb  of  the  sun  at 
the  Sacramento  Peak  Observatory  and  elsewhere . 

(b)  The  magnetic  field  survey  on  the  photospherlc  level  at  the  Kitt  Peak 
and  the  Mt .  Wilson  Observatories . 

(c)  If  available,  satellite  soft  x-ray  photographs  of  the  sun  (like 
Skylab  photographs) . 
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The  existing  network  of  solar  flare  patrol  systems  can  provide  the  basic 
input  data.  A  sudden  disappearance  of  chromospheric  filaments  appears  to  be 
another  (perhaps  more  Important)  source  (or  indication)  of  the  occurrence  of  a 


high  speed  stream  (Joselyn  and  McIntosh,  1981). 

(2)  Basic  codes 

(a)  Basic  solar  wind  code 

The  basic  solar  wind  code  (i)  analyzes  the  input  data  (1)  and  determines 
the  solar  magnetic  equator  on  a  sphere  of  2.5  solar  radii  (the  source 
surface);  (li)  infers  the  solar  wind  speed  distribution  on  the  source  surface; 
and  (iii)  computes  the  expected  solar  wind  speed  V,  density  n,  magnetic  field 
magnitude  B  and  orientation  (0,  d> )  at  the  location  of  the  earth. 

(b)  6  computation  code 

The  e  computation  code  takes  the  solar  wind  quantities  (V,  B,  9,  <t>)  at 
the  location  of  the  earth  as  a  function  of  time,  computed  by  the  basic  solar 
wind  code  and  computes  the  power  e(t)  generated  by  the  solar  wind- 
magnetosphere  generator  as  a  function  of  time. 

(c)  Auroral  oval  geometry  code 

The  oval  geometry  code  takes  the  power  e  as  the  input  and  determines  the 
geometry  of  the  poleward  boundary  of  the  auroral  oval  in  magnetic  latitude  - 
MLT  coordinates .  This  code  can  predict  also  the  equatorward  boundary  of  the 
oval,  along  which  the  auroral  E  appears. 


(d)  AS,  Dst  prediction  code 

This  code  takes  the  power  e  as  the  input  and  predicts  the  AE  and  Dst 
indices . 

(e)  Magnetospheric  convection  code 

This  code  takes  the  power  e  as  the  input  and  determines  the  convection 
pattern  of  magnetospheric  plasma  over  the  entire  polar  ionosphere.  It  has 
been  determined  empirically'  that  the  cross-polar  cap  poten'  ' al  is  proportional 
to  /e  (Reiff  et  al.,  1981).  Combining  this  empirical  re;  t  with  the  auroral 
oval  geometry  code,  one  can  infer  the  distribution  of  t'  alectrlc  field  in 
the  polar  regions. 

(f)  Empirical  auroral  oval-trough  code 

This  code  provides  the  geometry  of  the  oval  and  the  trough  for  given  AE 
and  Dst  values  as  the  input  data.  The  basis  for  this  code  is  an  empirical 
study  of  the  relationship  between  the  oval  geometry  and  the  AE,  Dst  indices. 
At  the  present  time,  the  oval  geometry  can  be  provided  graphically  for  a  given 
Q  index  or  AE  index.  This  code  is  a  backup  to  the  auroral  oval  geometry  code. 

(g)  Ionospheric  code 

This  code  determines  the  electron  density  profile  of  the  F2  region  as  a 
function  of  latitude  and  KLT  for  a  given  magnetospheric  convection  flow 
condition  on  a  given  day  by  using  the  magnetospheric  convection  code  as  the 
input  data  (Watkins,  1978;  Sojka  et  al.,  1982). 
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(h)  Energetic  particle  drift  code 

Since  the  mantle  aurora  is  produced  by  energetic  electrons  which  drift 


eastward  in  the  morning  sector,  this  code  is  needed  in  determining  the 
location  of  the  mantle  aurora, 

(i)  Prediction  of  the  equatorward  boundary  of  the  auroral  oval  and  trough 

The  codes  (c),  (f)  and  (g)  can  be  developed  in  such  a  way  that  their 

outputs,  namely  the  equatorward  boundary  of  the  oval  (of  both  discrete  and 
diffuse  auroras)  and  the  ionospheric  trough,  are  displayed  on  a  computer 
display  screen  on  a  continuous  basis. 

(j)  Flare-generated  wind  code 

The  scheme  presented  in  the  above  provides  the  'background'  conditions  of 
the  solar  wind  and  the  auroral  oval,  etc.  A  solar  flare  introduces  an 
additional  disturbance  to  the  background  condition.  The  solar  wind  condition 
after  a  solar  flare  depends  greatly  on  the  solar  wind  condition  prior  to  a 
flare  (namely,  the  background  condition)  since  the  former  results  partly  from 
the  interaction  between  the  slow  background  wind  and  the  flare-generated  fast 
wind . 

The  flare  wind  code  takes  (i)  the  observed  location  of  the  flare  with 
respect  to  the  solar  magnetic  equator,  (ii)  the  central  meridian  distance; 
(ill)  its  time  development;  and  (iv)  the  maximum  intensity  as  the  input 
data.  These  input  data  enable  us  to  introduce  "disturbance"  into  the 
background  basic  solar  wind  code,  simulating  numerically  disturbed  solar  wind 
conditions.  A  sudden  disappearance  of  chromospheric  filaments  appears  to  be 
another  (perhaps  more  important)  Indication  of  the  occurrence  of  a  high  speed 


stream.  However,  It  Is  necessary  to  quantify  the  observations  In  terms  of 
(111)  and  (lv). 


II.  Establishment  of  a  practical  geomagnetic  storm  prediction  scheme 


Summary 

He  have  contributed  in  establishing  a  direct  communication  link  between 
the  ISEE/C  satellite  and  the  Geophysical  Institute  via  the  Space  Environmental 
Laboratory,  Boulder,  Colorado  in  predicting  geomagnetic  storm  activity.  This 
direct  link  has  been  very  useful  for  the  operation  of  rocket  launching  at  the 
Poker  Flat  Rocket  Range.  In  fact,  the  system  was  operated  for  the  first  time 
in  1981  when  Air  Force  rockets  were  launched. 

During  the  last  several  years,  we  have  emphasized  the  Importance  of 
monitoring  this  solar  wind  parameter  on  a  real  time  basis  by  a  satellite  which 
Is  located  at  the  front  of  the  earth.  Fortunately,  the  International  Sun- 
Earth-Explorer  C  was  placed  at  the  llbratlon  point,  about  1.6  million 
kilometers  from  the  earth  toward  the  sun,  in  the  fall  of  1978.  From  its 
location,  the  satellite  monitors  the  solar  wind  and  relays  the  data  to  ground 
tracking  stations,  then  to  the  NASA  Goddard  Space  Flight  Center  (Maryland), 
and  to  the  NOAA  Space  Environmental  Laboratory  (Colorado)  where  the  solar 
wind-magnetosphere  energy  coupling  function  e  is  computed.  As  an  experiment, 
the  computed  data  were  transmitted  to  the  GEOS  satellite  over  the  Pacific 
Ocean  and  finally  relayed  to  the  Geophysical  Institute  during  the  spring 
rocket  campaign  In  1981.  It  takes  about  60  minutes  for  the  solar  wind  to 


reach  the  earth  from  the  satellite  location,  while  it  takes  only  about  7 
seconds  for  data  transmission.  Figure  2  shows  schematically  the  communication 
link  between  the  ISEE/C  satellite  and  the  Geophysical  Institute.  Figure 
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3  shows  our  receiving  Instrument.  Because  of  a  limited  communication 
capability,  e  values  are  coded;  2  x  1018  -  2,8,  5  x  lO1^  -  5,9,  7  x  1017  - 
7,7.  The  coded  values  were  sent  from  Boulder  every  5  minutes.  When  e  values 
were  less  than  10*8,  auroras  became  quiet.  When  e  values  exceeded  10^  for  an 
extended  period  (~3  hrs),  auroras  remained  In  the  southern  sky,  unsuitable  for 
rocket  launches.  When  e  values  increased  from  very  low  values  (<10*7  erg/sec) 
to  10  erg/sec  or  greater,  auroral  activity  increased  significantly  about  40 
minutes  afterwards.  Thus,  this  experiment  was  quite  successful  and  useful  for 
rocket  scientists,  since  they  were  able  to  schedule  the  launch  of  their 
rockets  in  time  to  intercept  the  forecasted  auroral  activity.  The  Global 
Weather  Central,  Offutt  AFB  is  also  monitoring  e  on  a  real  time  basis  by  using 
solar  wind  quantities  which  are  telemetered  from  the  ISEE/C  satellite  at  the 
llbration  point  (located  at  a  distance  of  240  earth  radii  toward  the  sun) . 

III.  Computer  display  system  of  the  auroral  oval 

We  set  up  a  computer  scheme  to  display  the  location  of  the  auroral  oval 
for  a  given  Q  at  a  given  UT  on  a  CRT  with  the  polar  geographic  map.  This  is  a 
computer  version  of  the  templet  scheme  constructed  by  Whalen  (1970)  "Auroral 
oval  plotter  and  monograph  for  determining  corrected  geomagnetic  local  time, 
latitude,  and  longitude  for  high  latitudes  in  the  northern  hemisphere." 
Figures  4a-d  show  some  of  the  examples  taken  from  the  CRT  for  AE  -  Oy,  140y, 
400y  and  660y  at  00  UT. 

We  have  also  constructed  a  model  of  the  magnetosphere  which  depends  on 
the  interplanetary  magnetic  field  (IMF).  It  is  possible  to  determine  the 
"open  region"  in  the  polar  cap  for  the  IMF  of  a  given  magnitude  and 
orientation.  Assuming  that  the  boundary  of  the  open  region  delineates 


approximately  the  auroral  oval,  we  have  devised  a  computer  scheme  to  show  the 


open  area  on  CRT,  An  Important  remaining  taek  la  to  convert  the  Q  Index  to 
the  AS  Index. 


W.  Prediction  of  the  development  of  geomagnetic  atorma  using  the  aolar  wlnd- 
magnetoaphere  energy  coupling  function  e 

Summary 

In  this  section,  we  explain  how  the  Dst  and  AS  prediction  codes  have  been 
developed.  In  particular,  It  la  shown  that  by  monitoring  time  variations  of 
the  solar  wlnd-magnetosphere  energy  coupling  function  e(t)  upstream  of  the 
solar  wind,  one  should  be  able  to  predict  fairly  accurately  the  growth  and 
decay  of  individual  magnetospherlc  substorms  and  storms.  In  this  section,  we 
show  how  the  two  key  geomagnetic  atom  Indices,  Dst  and  AS,  can  be  estimated 
on  the  basis  of  the  predicted  values  of  e.  For  details  of  this  study,  see 
Akasofu  (1981,  1982) 

1 .  Formulation 

It  is  generally  agreed  that  the  major  part  of  the  energy  transferred  from 
the  solar  wind  to  the  magnetosphere  is  deposited  In  the  ring  current,  so  that 
we  assume  that  70%  of  e  is  consumed  in  the  formation  of  the  ring  curreat  belt, 
namely  that  the  rate  UR  of  ring  current  particle  Injection  Is  given  by 
UR  *  0.7e.  The  predicted  ring  current  field  DR(y)  is  then  obtained  by 
integrating  the  equation 

.  ,_20  ,508  .  DRn  a  _ 

4  x  10  (-r-  +  — )  -  0.7e 

St  tr 

where  xR  Is  the  lifetime  of  the  ring  current  particles .  The  dependence  of  x., 


on  storm  intensity  is  taken  to  be: 


■  20  h  for  e  <  10*®  ergs  1 
-  2h  for  1018  ergs'1  e  <  1019  ergs'1 

■  1  h  for  e  >  1019  ergs'1 

As  we  shall  see  later,  this  set  of  values  of  xR  tends  to  predict  slightly 
greater  aain  phase  decreases  than  observed  ones.  However,  we  shall  use  this 
set  in  this  paper,  since  a  similar  set  has  been  used  extensively  in  the 
past . 

The  Dst  index  consists  of  at  least  two  components,  the  ring  current  field 
DR(y)  and  the  solar  wind  compression  field  DCF(y) .  The  solar  wind  compression 
field  is  given  by 

DCF(y)  •  13.1  x  104(/^  - 

where  Pj  and  P2  denote  the  solar  wind  pressure  before  and  after  storm  onset. 
Therefore,  taking  into  account  (2)  and  (3),  the  predicted  Dst  index  (Dstp)  is 
given  by 

Dstp  -  DCF  -  DR. 


In  the  next  section,  the  baseline  of  Dstp  is  adjusted  to  give  Dst  -  Dstp  at 
the  beginning  of  each  storm  period. 


There  Is  still  much  uncertainty  between  the  Joule  heat  production  rate  Uj 
In  the  Ionosphere  and  the  AE  index,  although  the  rate  Uj  in  the  polar 
ionosphere  is  10-30Z  of  the  power  e  generated  by  the  solar  wind-magnetosphere 
dynamo.  One  of  the  complications  Involved  in  this  particular  study  is  that  Uj 
is  approximately  proportional  to  e  for  small  values  of  e  (<10  ergs  ),  but 
this  relation  breaks  down  for  e  >  10  ergs  .  Therefore,  we  use  here  an 
empirical  relationship  between  e  and  the  AE  index  which  includes  the  fact  that 

1  Q  _1 

the  AE  index  tends  to  "saturate”  for  e  >  10  ergs  .  The  empirical 
expression  is  given  by 


AEp(Y)  -  -300  (log)2  +  11,700  log  c  -  11,3200. 


2.  Examples 

In  this  section,  we  study  first  one  of  the  simplest  storm  developments  in 
detail  and  then  examine  how  well  the  solar  wind-magnetosphere  coupling 
function  e  can  predict  the  development  of  geomagnetic  storms  in  terms  of  the 
predicted  Dst  index  (Dstp) . 

Figure  $  shows,  from  the  top,  the  solar  wind-magnetosphere  energy 
coupling  function  e,  Dstp,  the  observed  Dst,  AEp  and  the  observed  AE  index, 
for  the  16-18  March,  1974  storm.  No  time  delay  is  assumed  in  the  energy 
coupling  and  dissipation  processes. 

Prior  to  the  onset  of  the  storm  of  March  16,  1974  (sac  at  1128  UT),  there 

lO 

occurred  minor  substorm  activity  with  e  ~  2  x  10  ergs  ,  which  caused  a  weak 
ring  current,  as  indicated  by  a  slight  depression  of  the  Dst  index.  This 
change  is  reproduced  fairly  well  in  the  predicted  Dst  (Dstp).  At  storm  onset 
(1128  UT),  a  large  increase  of  the  solar  wind  pressure  caused  the  Dst  index  to 


Increase.  The  same  trend  can  be  seen  in  Dst  ,  although  it  did  not  become 


positive;  this  difference  Is  due  to  the  fact  that  the  intensity  of  the 
predicted  ring  current  field  was  slightly  greater  than  the  observed  one. 


The  ssc  was  immediately  followed  by  a  rectangular  wave-like  increase  of 
€,  lasting  for  about  8  h.  One  can  easily  see  that  the  predicted  Dst  variation 
(Ostp)  is  quite  similar  to  the  observed  one.  in  terms  of  not  only  the 
magnitude  of  the  decrease,  but  also  the  time  development.  The  ring  current 
field  developed  rapidly  until  about  the  time  when  the  rectangular  change 
subsided  and  then  decayed  slowly.  A  slight  irregularity  in  the  decay  was 
associated  with  minor  substorm  activity  on  17  and  18  March.  The  AE  index 
variation  during  the  storm  consisted  of  a  main  pulse  associated  with  the 
rectangular  wave-like  change  of  e  and  smaller  Increases  before  and  after  it . 
The  main  pulse  was  reasonably  well  reproduced,  but  the  two  smaller  changes 
were  exaggerated  a  little  in  the  predicted  AE  index  (AEp) . 

2.1  Storm  with  an  Insignificant  main  phase 

In  this  type  of  storm,  a  significant  main  phase  falls  to  develop,  in 
spite  of  a  distinct  ssc  and  the  initial  phase.  The  Dst  index  in  Figure  6 
shows  that  the  storm  of  13-15  February  1973  is  of  this  type.  It  is  easy  to 
see  that  this  particular  solar  wind  flow  contained  only  a  small  amount  of  e. 
The  predicted  Dst  variation  is  quite  similar  to  the  observed  one.  If  the 
baseline  can  be  adjusted  for  an  increase  of  the  observed  Dst  at  about  20  UT  on 
13  February,  the  agreement  between  Dst  and  Dstp  can  become  even  better. 

2.2  Mini-storm 

This  type  of  storm  develops  as  a  result  of  a  very  Intense  Impact  of  the 
storm  (causing  an  unusually  large  ssc),  but  the  main  phase  is  weak  and  its 
duration  is  short  (against  the  expectation  from  the  magnitude  of  ssc).  Figure 
7  shows  an  example  of  such  storms.  The  Intense  impact  of  the  solar  wind  can 
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clearly  be  recognized  by  a  large  ssc.  One  can  easily  recognize  that 

complicated  time  variations  of  the  observed  Dst  are  fairly  well  reproduced  by 

Dst„.  The  difference  between  the  observed  Dst  and  Dst  prior  to  storm  onset 
P  P 

Is  due  to  the  fact  that  the  substorm  activity  from  ~  22  UT  5  March  to  06  UT  6 
March  caused  an  appreciable  depression  of  Dstp,  but  this  depression  was  not 
registered  In  the  observed  Dst . 

2 .3  Typical  storms 

In  Figures  8  and  9  we  examine  typical  storms .  The  storm  of  7-8  February 

1967  began  with  a  rapid  Increase  of  e  at  the  time  of  ssc.  The  main  phase 

developed  rapidly  for  about  10  h  and  began  to  decay  near  the  end  of  7  February 

as  e  decreased.  The  Increase  of  e  at  about  06  UT  on  8  February  caused  a  new 

Intensification  of  the  main  phase.  One  can  clearly  see  that  the  development 

and  decay  of  the  main  phase  decrease  Is  fairly  well  reproduced  by  the 

predicted  Dst .  The  AE  Index  Is  also  fairly  well  reproduced . 

In  Figure  9,  the  storm  of  31  March  -  2  April  1973  began  with  a  medium 

increase  of  e  at  about  15  UT  on  31  March,  which  caused  a  weak  main  phase. 

This  increase  of  e  ended  at  about  22  UT  on  the  same  day.  Another  Increase  of 

19  —1 

e  ended  at  about  23  UT,  but  its  magnitude  was  less  than  10  ergs  until 

about  12  UT  on  1  April.  As  a  result,  the  Dst  decrease  remained  almost 

constant  for  more  than  12  h.  A  new  increase  of  e  began  at  about  09  UT  on  1 

19 

April.  An  Intense  ring  current  belt  began  to  develop  when  e  reached  ~  10 
ergs-* .  The  main  phase  began  to  decay  near  the  end  of  1  April.  Again,  the 
overall  pattern  of  the  development  of  the  main  phase  decrease  is  well 
reproduced  by  Dstp.  The  predicted  AE  index  agrees  also  fairly  well  with  the 
observed  one.  However,  during  the  period  when  the  main  phase  was  rapidly 
developing  the  observed  AE  Index  decreased . 
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Figure  10  shows  an  example  In  which  the  predicted  Dst  was  much  greater 
than  the  observed  Dst,  although  the  overall  feature  of  the  time  development 
was  fairly  well  reproduced.  Jn  this  storm  of  15-18  ^ebrufry,  1967,  e  reached 
one  of  the  highest  values  during  the  last  15  years.  There  are  at  least  two 
reasons  for  this  discrepancy  between  the  observed  and  predicted  Dst 
variations.  The  first  is  that  tg  would  have  become  less  than  1  h  for  e  >  5  x 
1019  ergs-*  or  that  e  would  have  been  overestimated  by  assuming  that  JtQ  is  a 
constant.  A  large  increase  of  the  solar  wind  pressure  could  have 
significantly  reduced  AQ. 

2.4  Long-lasting  storms 

Medium  intensity  storms  last  for  about  48  h  or  so  and  decay  slowly.  Such 
storms  are  associated  with  a  long-lasting  increase  of  e,  with  rapid  changes, 

|Q  _| 

although  the  maximum  magnitudes  of  e  remain  less  than  10  ergs  .  Figure  11 
shows  an  example  of  such  a  storm.  Unlike  some  of  the  major  storms  which 
develop  rapidly  and  decay  rapidly  as  well,  this  particular  storm  was 
characterized  by  a  long-lasting  Dst  decrease,  although  the  magnitude  of  the 
Dst  values  was  not  very  large.  Again,  the  predicted  Dst  reproduces  fairly 
well  the  overall  features  of  this  type  of  storm. 

3.  Conclusion 

He  have  demonstrated  that  one  can  make  a  fairly  accurate  prediction  of 
the  time  development  of  the  main  phase  decrease  as  manifested  in  the  two 
geomagnetic  Indices  Dst  and  AE,  by  monitoring  e(t)  upstream  of  the  solar 
wind.  Satellite  ISEE/C,  anchored  at  the  llbratlon  point  (at  a  distance  of 
about  1.6  x  10^  km  toward  upstream  of  the  solar  wind),  can  provide  us  with  the 
necessary  information  to  allow  for  a  1  h  lead-time  forecasting.  A  longer  lead 
time  forecasting  would  require  numerical  forecasting  of  e  in  a  way  similar  to 
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numerical  weather  forecasting.  For  this  purpose,  one  must  be  able  to  simulate 
solar  wind  disturbances  to  determine  V  and  B  as  a  function  of  location  and 
time  in  Interplanetary  space,  and  computed  e  at  the  location  of  the  earth. 

V.  Dependence  of  the  geometry  of  the  region  of  open  field  lines  on  the 
interplanetary  magnetic  field 


Summary 

In  this  section,  we  explain  how  the  auroral  oval  geometry  code  has  been 
developed.  In  particular,  the  geometry  of  the  open  flux  area  in  the  polar 
region  is  computed  by  superposing  a  uniform  Interplanetary  magnetic  field 
(IMF)  with  various  orientation  angles  to  a  model  of  the  magnetosphere.  For 
details  of  this  study,  see  Akasofu  et  al .  (1981).  The  results  show  that  the 
IMF  By  component  is  as  Important  as  the  Bz  component  in  "opening"  the 
magnetosphere .  It  is  also  shown  that  the  computed  area  of  open  field  lines  is 
remarkably  similar  to  the  observed  ones  which  were  determined  by  using  the 
entry  of  solar  electrons.  In  particular,  when  the  IMF  vector  is  confined  In 
the  X-Z  plane  and  the  Bz  component  has  a  large  positive  value,  the  open  area 
becomes  crescent-shaped,  coinciding  approximately  with  the  cusp  region. 

1.  Introduction 

There  are  now  several  observations  which  indicate  or  suggest  that  drastic 
changes  of  the  geometry  of  the  open  region  occur  as  the  orientation  of  the 
Interplanetary  magnetic  field  (IMF)  vector  changes.  These  observations  are 
based  on  the  geometry  of  the  area  into  which  low-energy  or  high-energy  solar 
electrons  enter.  Here,  we  examine  specifically  the  dependence  of  the  geometry 
of  the  open  field  line  region  of  the  magnetosphere  on  the  IMF  by  assuming  a 
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simple  linear  superposition  of  Che  IMP  on  a  model  magnet os pheric  field.  This 
study  can  also  be  considered  Co  be  a  qualitative  test  of  the  validity  of  the 
superposition  method,  since  one  can  compare  the  observed  open  region  and  the 
computed  open  region.  It  is  indeed  of  great  interest  to  examine  whether  or 
not  this  simple  method  can  provide,  as  a  first  approximation,  a  topologically 
correct  configuration  of  the  open  magnetosphere,  without  using  a  full-scale 
plasma  simulation  method. 

In  order  to  make  the  open  region  of  the  model  as  realistic  as  possible, 
one  must  carefully  construct  a  three-dimensional  model  of  the  magnetosphere  so 
as  to  be  able  to  determine  whether  or  not  individual  magnetospherlc  field 

lines  are  connected  to  the  IMF  field  lines.  Our  model  of  the  magnetosphere 
consists  of  the  earth's  dipole  field  (of  magnetic  moment  Mg)  and  an  image 
dipole  field  corresponding  to  a  dipole  of  moment  7.0  Mg  at  X  -  28.0  Rg,  plus 
the  magnetic  field  of  the  magnetotail  (solar-magnetospherlc  coordinates  are 
used  in  this  paper) .  The  tail  field  is  simulated  with  two  sets  of  77  ring 
currents,  one  set  in  each  hemisphere,  separated  by  0.2  Rg  across  the 

equatorial  plane  (the  so-called  "neutral  sheet").  The  currents  flow  counter¬ 
clockwise  (as  viewed  from  the  earth)  in  the  northern  ring  and  clockwise  in  the 
southern  ring.  They  are  distributed  uniformly  over  a  distance  of  180  Rg,  the 
nearest  pair  being  located  at  X  ■  -10  Rg;  thus  the  successive  rings  are  spaced 
uniformly  by  a  distance  of  2.37  Rg.  In  order  to  generate  the  observed  8° 

flaring  of  the  tall,  the  radius  of  the  rings  is  made  to  Increase  linearly  by 

14Z.  Further,  the  current  Intensity  in  successive  rings  is  made  to  decrease 
with  distance  down  the  tail,  starting  with  250,000  A  at  X  -  -10.0  Rg  to  zero 
at  X  -  -190  Rg. 

In  this  way,  we  construct  first  a  model  of  the  magnetosphere  in  the 
absence  of  the  IMF  and  determine  the  geometry  of  the  magnetopause.  The 


magnetopause  consists  of  the  magnetic  field  lines  which  "pass**  through  a  very 
small  area  around  the  dayside  neutral  points  (one  in  each  hemisphere).  Then, 


superposing  the  IMF  of  magnitude  |B|  -  10y  in  the  model,  we  trace 

systematically  a  large  number  of  field  lines  in  the  magnetosphere.  Following 

Saunders  (1976),  the  open  field  lines  are  defined  as  magnetic  field  lines 

which  reach  the  magnetopause.  The  use  of  a  relatively  large  value  of 


|  B |  ■  10y  as  the  magnitude  of  the  IMF  is  to  examine  effects  of  the  IMF  for  the 
magnetosphere  with  a  relatively  short  tall  of  length  200  Rg.  It  is  generally 
accepted  that  the  length  of  the  tall  is  significantly  longer  than  200  Rg,  so 
that  a  smaller  magnitude  of  the  IMF  will  have  the  same  effects . 


2.  Results  and  comparison  with  observations 


The  geometry  of  the  open  field  line  area  is  determined  for  two  sets  of 
IMF  orientations.  In  the  first  set,  the  IMF  vector  is  confined  In  the  X-Z 
plane,  and  the  angle  between  the  Z-axis  and  the  IMF  vector  is  denoted  by  a, 
being  reckoned  positive  as  the  vector  rotates  away  from  the  Sun.  The  values 
of  a  examined  here  are  30°,  60°,  90°,  120°  and  180°,  so  that  the  Bx  component 
is  always  negative.  In  the  second  set,  the  IMF  vector  is  confined  in  the  T-Z 
plane,  and  the  angle  between  the  Z-axls  and  the  IMF  vcector  is  denoted  by  8, 
being  reckoned  positive  as  the  vector  rotates  toward  the  dawn  side,  so  that 
the  By  component  is  always  negative.  The  values  of  8  examined  here  are  also 
30°,  60°,  90°,  120°  and  180#. 

(a)  a  -  30°  (Bx  -  -5.0y,  By  -  Oy,  Bz  -  8.7y) 

P  “  30°  (Bx  -  0y,  By  -  -5.0y,  Bz  -  8.7y). 

Figure  12  shows  the  geometry  of  the  northern  open  region  for  both  sets  of 
IMF  orientations,  a  -  30°  and  8  ■  30°  in  latitude-magnetic  local  time  (MLT) 
coordinates.  For  a  •  30°,  the  open  region  is  a  very  small  area  which  is 


located  at  79°  in  latitude  along  the  noon  meridian.  Therefore,  the  magneto¬ 
sphere  is  practically  closed  in  this  situation. 

For  P  »  30°,  the  open  region  is  almost  confined  in  the  evening  sector  in 
the  northern  hemisphere.  Though  it  is  not  shown  here,  the  open  region  is 
located  in  the  morning  sector  in  the  southern  hemisphere  and  its  geometry  is 

the  mirror  image  (with  respect  to  the  noon-midnight  meridian)  of  the  northern 

one.  The  significance  of  this  result  will  he  discussed  in  (b) . 

(b)  a  -  60°  (Bx  -  -8.7y,  By  -  0y,  \  -  5.0Y) 

P  -  60°  (Bx  -  0y,  By  -  -8.7y,  Bz  -  5.0y). 

Figure  13  shows  the  geometry  of  the  northern  open  region  for  both  A  **  60° 
and  B  »  60°.  One  can  immediately  recognize  a  large  difference  of  the  geometry 
of  the  open  region  for  the  same  value  (60°)  of  a  and  p. 

For  a  -  60°,  the  magnetic  field  lines  from  a  crescent-shaped  region  in 
the  dayside  can  reach  the  magnetopause.  It  is  interesting  to  note  in  this 
respect  that  HcDlarmld  et  al.  (1980)  found  that  high-energy  solar  electron 
fluxes  have  a  single  peak  in  their  pitch-angle  distribution  in  a  crescent¬ 
shaped  area  in  the  polar  region  and  a  double  peak  in  the  rest  of  the  polar 

region  when  the  IMF  components  were  B„  ~  -20y,  B_ - 5y  and  B_  ~  +30y.  They 

suggested  that  the  open  region  had  a  crescent  shape.  Their  figures  are 
reproduced  here  as  Figures  14(a)  and  (b).  One  can  see  the  similarity  of  the 
open  regions  in  Figure  14(a)  and  Figure  13  (a  -  60°).  It  Is  interesting  to 
note  also  in  this  connection  that  the  dayside  auroral  arc  system  lies  in  a 
crescent  region  similar  to  those  shown  in  Figures  13  and  14(a)  and  remains 
bright  even  when  the  IMF  Bz  component  is  large  and  positive  (although  the 
nlghtslde  arc  system  becomes  very  faint). 

For  P  -  60°,  the  open  region  is  mainly  confined  to  the  evening  sector  in 
the  northern  hemisphere;  in  the  southern  hemisphere,  the  open  region  is  mainly 


confined  Co  Che  morning  secCor.  Again,  Che  corapuCed  open  region  la  similar  Co 

ChaC  suggesced  by  McDiarmid  eC  al.  (1980),  alchough  Che  Bx  component  was 

greaCer  than  the  By  component,  namely  Bx  =  4.3y,  By  ■  -0.7y  and  Bz  =■  +9.2y. 

The  open  area  in  the  southern  hemisphere  Is  Che  mirror  Image  (with  respect  to 

the  noon-midnight  meridian)  of  it  and  appears  In  the  morning  sector.  If  the 

IMF  vector  is  directed  In  the  opposite  direction,  namely  B  ■  0y,  B  -  8.7y 

x  y 

and  Bz  =■  5.0y.  Such  an  asymmetry  of  the  entry  of  solar  electrons  with  respect 
to  the  sign  of  the  IMF  By  component  was  noted  by  Meng  et  al,  (1977),  Further, 
Meng  et  al.  (1977)  showed  that  the  situation  is  reversed  in  the  southern 
hemisphere.  Therefore,  our  results  are  qualitatively  in  agreement  with  the 
observations . 

(c)  a  -  90°  (Bx  -  -10y,  By  -  0y,  Bz  -  0y) 

P  =  90°  (Bx  -  Oy,  By  -  -10y,  Bz  -  0y) . 

Figure  15  shows  the  results  for  a  =■  90°  and  p  ■  90°.  The  results  for  a  * 
90°  are  similar  to  those  for  a  -  60°,  except  that  the  open  areas  become 
appreciably  greater  than  those  for  a  *  60®.  It  is  thus  of  great  interest  to 
examine  auroras  under  similar  IMF  conditions.  Therefore,  the  nightside  of  the 
polar  region  is  still  essentially  closed  even  for  a  *  90°.  For  P  *  90°, 
however,  the  open  area  expands  drastically  and  becomes  similar  to  that 
enclosed  by  the  average  size  auroral  oval.  This  result  indicates  that  the  By 
component  is  very  important  in  making  the  nightside  of  the  polar  region 
open.  For  IMF  polar  angles  of  less  than  90®  or  B_  >  0,  the  importance  of  the 
iy  component  in  opening  the  magnetosphere  is  indeed  very  crucial. 

(d)  a  -  120®  (Bx  -  -8.7y,  By  -  0y,  Bz  -  -5.0y) 

P  -  120®  (Bx  -  0y,  By  -  -8.7y,  Bz  -  -5.0y) 

Figure  16  shows  results  for  a  »  120®  and  p  -  120®.  For  both  cases,  the 
open  area  Is  close  to  the  region  enclosed  by  the  normal  size  oval.  This  is 
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expected  because  the  IMF  B_  component  is  negative.  Note  that  there  la  a 
slight  asymmetry  of  the  area  with  respect  to  the  noon-midnight  meridian  for 


P  -  120°. 

(e)  a  -  P  -  180°  (Bx  -  0y,  By  -  0y,  Bz  -  -10y). 

Figure  17  shows  the  result.  The  auroral  oval  is  the  same  for  both  a  - 
180°  and  p  *  180°.  The  open  area  is  a  little  greater  than  that  of  the  normal 
size  oval. 

3.  Discussion 

The  remarkable  similarity  of  the  open  regions  (inferred  from  the  Impact 
area  of  solar  electrons)  by  McDlarmid  et  al.  (1980)  and  Meng  et  al.  (1977)  and 
the  computed  open  regions  proves  the  validity  of  the  simple  superposition 
method  of  the  IMF  on  a  magnetospherlc  field  model,  at  least  in  studying  the 
topological  changes  of  the  magnetospherlc  structure  caused  by  variations  of 
the  orientation  of  the  IMF. 

VI.  Polar  cap  arcs  over  Thule 

The  need  for  a  computerized  analysis  of  all-sky  auroral  photographs  has 
been  stressed  during  the  last  several  years.  Therefore,  we  have  established  a 
scheme  to  analyze  all-sky  photographs  by  using  a  digitizer  and  to  display  the 
location  of  auroral  arcs  on  a  geographic  (or  geomagnetic)  map  and  in  Invariant 
latitude  -  MLT  coordinates.  The  display  is  accomplished  by  specifying  the 
geographic  location  of  a  fixed  or  moving  observatory  and  the  universal  time 
(UT)  when  the  photographs  are  taken  as  well  as  the  type  of  camera  used.  The 
lens  characteristics  are  digitized  for  the  Fairchild  lens,  the  University  of 
Alaska  all-sky  (fish  eye)  camera  and  the  IGY  (mirror  system)  all-sky  camera. 


The  scheme  has  been  applied  to  all-sky  photographs  taken  during  the  1981 
Thule  expedition  by  the  USAF  Flying  Observatory.  Figure  18a  shows  some  of  the 
auroral  photographs  which  were  taken  on  January  21,  1981,  showing  a  rapid 
development  of  a  polar  cap  arc  over  Thule.  Figure  18b  shows  the  location  of 
the  arc  In  Figure  18a  In  Invariant  latitude  -  MLT  coordinates . 
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